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We determine charge transfer statistics in a quantum conductor driven by a time-dependent 
voltage and identify the elementary transport processes. At zero temperature unidirectional and 
bidirectional single charge transfers occur. The unidirectional processes involve electrons injected 
from the source terminal due to excess dc bias voltage. The bidirectional processes involve electron- 
hole pairs created by time-dependent voltage bias. This interpretation is further supported by 
the charge transfer statistics in a multiterminal beam splitter geometry in which injected electrons 
and holes can be partitioned into different outgoing terminals. The probabilities of elementary 
processes can be probed by noise measurements: the unidirectional processes set the dc noise level 
while bidirectional ones give rise to the excess noise. For ac voltage drive, the noise oscillates 
with increasing the driving amplitude. The decomposition of the noise into the contributions of 
elementary processes identifies the origin of these oscillations: the number of electron-hole pairs 
generated per cycle increases with increasing the amplitude. The decomposition of the noise into 
elementary processes is studied for different time-dependent voltages. The method we use is also 
suitable for systematic calculation of higher-order current correlators at finite temperature. We 
obtain current noise power and the third cumulant in the presence of time-dependent voltage drive. 
The charge transfer statistics at finite temperature can be interpreted in terms of multiple charge 
transfers with probabilities which depend on energy and temperature. 



PACS numbers: 72.70.+m, 72.10.Bg, 73.23.-b, 05.40.-a 



I. INTRODUCTION 



The charge transmitted through a mesoscopic conduc- 
tor during a fixed time interval fluctuates because of 
charge discreteness and stochastic nature of transport. 
The objective of the statistical theory of quantum trans- 
port, full counting statistics, is to completely characterize 
the probability distribution of transferred charge. The 
field of full counting statistics has attracted significant 
attention because it provides the most detailed informa- 
tion on charge transfer accessible in the measurements 
of the average current, current noise power and higher- 
order current correlations .^^i^ii^ 

The problem of full counting statistics has been ad- 
dressed first by Levitov and Lesovik^'*' for the case 
of dc biased multiterminal junctions and subsequently 
generalized to a time-dependent voltage biasiii^ In ad- 
dition to the scattering theory of Levitov and Leso- 
vik (see also Ref. 0), the theoretical approaches to full 
counting statistics include the so-called stochastic path- 
integral approacbiSiii and the quantum-mechanical the- 
ory based on an extension of the Keldysh-Green's func- 
tion technique ji^ii^iii^'i^ii^ Although equivalent, these 
theories provide different methods to access the charge 
transfer statistics. Using the extended Keldysh-Green's 
functions technique, the theory of full counting statis- 
tics of a general quantum-mechanical observable has 
been formulatedjiiii^ including the effects of a detec- 
tor back action.— The problem of quantum measure- 
ment in the context of full counting statistics of non- 
commuting spin observables has also been analyzedi^>^ 



The discretization^^'^'^ in space of the Keldysh-Green's 
functions results in a quantum circuit theory^ ^'^^i^^ 
which greatly simplifies the calculations. The circuit the- 
ory can describe junctions with different connectors and 
leads, as well as multiterminal circuits. 

The goal of evaluation of full counting statistics 
in mesoscopic conductors has essentially been accom- 
plished: the aforementioned theories provide methods 
and techniques to calculate the cumulant generating 
function of transferred charge from which the probabil- 
ity distribution can be obtained. However, the interpre- 
tation of the resulting statistics is often not straightfor- 
ward. This is because the total statistics are composed of 
many electrons injected towards conductor, each exhibit- 
ing different chaotic scattering and multiple reflections 
before entering an outgoing terminal. At finite temper- 
atures, the charge flow is bidirectional with both injec- 
tion and absorption of electrons from terminals. The 
particle correlations due to Pauli principle, interactions, 
and those induced by the superconducting and/or ferro- 
magnetic proximity effect also affect the charge transfer 
statistics. 

To understand the statistical properties of collective 
charge transfer, one has to find the independent ele- 
mentary processes constituting the statistics. The in- 
dependent processes can be revealed by decomposition 
of the total cumulant generating function into a sum of 
simpler ones. The physical interpretation obtained this 
way goes beyond information contained in the average 
current, noise, and finite-order cumulants, and pertains 
to all transport measurements. For example, the ele- 
mentary processes in a dc biased normal junction are 



2 



single-electron transfers, which are independent at dif- 
ferent energies. At finite temperatures, the electrons are 
transferred in both directions with correlated left- and 
right-transfers. At low temperatures only unidirectional 
electron transfers remain, with direction set by polar- 
ity of applied voltage. The elementary processes in a 
superconductor/normal contact are single- and double- 
charge transferSi^i^i At low temperatures and voltages 
below the superconducting gap only double-charge An- 
dreev processes remain, while above the gap normal 
single-charge transfers dominate. The elementary charge 
transfer processes between superconductors with dc bias 
applied consist of multiple charge transfers due to mul- 
tiple Andreev reflect ions i^^ii^'^i^ The interpretation of 
elementary charge transfer processes between supercon- 
ductors with constant phase difference is more subtle. 
In this case charge transfers acquire formally negative 
probabilities and the proper physical interpretation has 
to include dynamics of a detector. ^'^'—1^2, 

In contrast to dc voltage bias, time-dependent volt- 
age drive mixes electron states of different energies which 
in combination with the Pauli principle leads to a non- 
trivial charge transfer statistics. The signatures of this 
statistics have been studied first through the noise in 
an ac-driven junction (photon-assisted noise) i^^*^ The 
photon-assisted noise at low temperatures is a piecewise 
linear function of a dc voltage offset with kinks corre- 
sponding to integer multiples of the driving frequency 
and slopes which depend on the amplitude and shape of 
the ac component. This dependence has been observed 
experimentally in normal coherent conductor a^^'^^ and 
diffusive normal metal/superconductor junctionSi^ The 
photon-assisted noise at small driving amplitudes is due 
to an electron-hole pair which is created by ac drive (with 
a low probability per voltage cycle) and injected towards 
the scatterer.— i^i^ However, to obtain the elementary 
processes, the knowledge of full charge transfer statistics 
is needed. 

The elementary processes in the presence of time- 
dependent drive have been identified in Refs . [ilHi^lisI for 
the special choice of a driving voltage. The authors have 
studied Lorentzian voltage pulses of the same sign car- 
rying integer numbers of charge quanta. In this case the 
elementary processes are electrons injected towards the 
scatterer without creation of the electron- hole pairs. The 
charge transfer is unidirectional with binomial statistics 
set by the effective dc voltage. This result is independent 
of the relative position of the pulses, their duration, and 
overlap. The many-body quantum state generated by 
these pulses has been obtained recently. An alternative 
way to inject single electrons free of electron-hole pairs 
is by using time-dependent shifts of a resonant level in a 
quantum dot4i.i^.i&iil.ia 

The elementary charge transfer processes for arbi- 
trary time-dependent voltage drive have been identified 
in Ref. At low temperatures these processes are 

single-charge transfers which originate from electrons and 
electron-hole pairs injected towards the scatterer. The 




FIG. 1: The elementary processes shown in different regimes. 
For a dc bias voltage V and at arbitrary temperature Te, the 
electrons are transferred independently at different energies. 
The probability of transfer is given by occupation numbers 
fi;2{£) of the leads and transmissions {Tp} of the junction. 
The transport occurs when the incoming electronic state is 
occupied, the outgoing state is empty, and the electron is 
transferred across the junction (a). On the other hand, an ac 
voltage drive AV{t) mixes the electron states of different en- 
ergies. At zero temperature it creates the electron-hole pairs 
with probabilities pk (b), in addition to electrons injected due 
to excess dc voltage. A created electron-hole pair contributes 
to transport only if one particle from a pair is transmitted and 
the other is reflected, which occurs with probability TpRp {Rp 
being reflection coefficient). Finally, at finite temperature and 
in the presence of a time-dependent voltage drive, the statis- 
tics can be interpreted in terms of multiple n-charge transfers 
with probabilities q„ (c). 

electrons are injected due to excess dc voltage applied 
and give a binomial contribution to the total cumulant 
generating function. The electron-hole pairs are created 
by the ac component of the voltage drive. The proba- 
bilities of pair creations per voltage cycle depend on the 
details of the ac drive. For the special choice of opti- 
mal Lorentzian voltage pulses no electron-hole pairs are 
created. In general, however, an ac drive does create 
electron-hole pairs, with more and more pairs created per 
voltage cycle as the amplitude of the drive increases. A 
geometric interpretation of elementary processes is stud- 
ied in Ref. [sol . For time- dependent scatterer, the con- 
straints imposed on charge transfer statistics restrict the 
allowed charge transfersi^ 

In this article we present a comprehensive study of el- 
ementary processes in a generic mesoscopic conductor in 
different physical regimes, as depicted in Fig. [TJ We first 
obtain the cumulant generating function for a coherent 
voltage-driven quantum contact at finite temperatures, 
using both the formalism of Levitov and Lesovik and the 
circuit theory of quantum transport (Sec. |TT|. We then 
identify the elementary charge transfer processes at zero 
temperature fSec. lIIl"A| . Decomposition of the noise into 
contributions of elementary processes for different voltage 
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signals is studied in Sec. IIIIBI For ac voltage drive, the 
differential noise oscillates as the amplitude of the drive 
increases due to new electron-hole pairs being created 
per voltage cycle. The number of electron-hole pairs per 
cycle becomes large for driving amplitudes much larger 
than frequency. In this case the statistics in the leading 
order reduces to uncorrelated electron and hole transfers 
and depends on an effective voltage only, independent of 
the details of the ac drive. The elementary processes can 
be probed also for both ac and dc voltages present, e.g., 
in the regime in which the ac component of the drive is 
kept fixed and the dc offset changes. In this case both 
the electrons and created electron-hole pairs are injected 
towards the junction. The electron-hole pairs give rise to 
the excess noise with respect to the noise level set by dc 
voltage offset. 

The notion of elementary processes, being a mat- 
ter of interpretation of the full charge transfer statis- 
tics, provides physically plausible and intuitive descrip- 
tion of transport in terms of statistically independent 
events. Such a description persists in all transport mea- 
surements. For example, the picture of electrons and 
electron-hole pairs injected towards the scatterer, which 
we have inferred from a study of a 2-terminal junction, is 
further supported by the evaluation of the full counting 
statistics in a beam-splitter geometry in which case the 
injected particles can be separated into different outgo- 
ing terminals. These processes can be probed by current 
cross correlations fSec. IlllUj) . 

The method we use is also suitable for systematic cal- 
culation of the higher-order current correlators at finite 
temperature. The current noise power and the third cu- 
mulant in the presence of time-dependent voltage drive 
are obtained in Sec. IIVI The interpretation of the full 
counting statistics and the elementary processes are fun- 
damentally different at finite temperature. For a periodic 
drive, the elementary processes are transfers of multiple 
integer charge quanta with probabilities which depend on 
energy and temperature fSec. W\. 



II. CUMULANT GENERATING FUNCTION 

The system we consider is a generic 2-terminal meso- 
scopic conductor characterized by a set of transmission 
eigenvalues {Tp}. For a large number of channels, the 
transport properties are universal and independent of 
microscopic details of geometry of the junction and po- 
sitions of impurities. We can neglect the energy depen- 
dence of transmission eigenvalues if the the electron dwell 
time is small with respect to time scales set by the in- 
verse temperature and applied voltage. The cumulant 
generating function of the transferred charge is given by 
determinant formul a^i'^^i^^ 



where 



/i 

/2 



s = 



r t' 
t r' 



1 



(2) 



Here f is the matrix of occupation numbers of the termi- 
nals which is diagonal in the terminal indices and scalar 
in transport channels, S is the scattering matrix of the 
junction, and . . . A^^, is the transformation which in- 
corporates the counting fields. Because the current is 
conserved, it is sufficient to assign only one counting field 
X to the left terminal. We count charges which enter the 
left terminal irrespective the channel, the energy, or the 
spin, with x being scalar in the corresponding indices. 
The determinant in Eq. ([T]) is taken with respect to the 
terminal, channel, energy, and spin indices. 

The occupation numbers /i and /2 of the left and 
right terminals arc matrices in energy indices and scalars 
in channel and spin indices. In the case of dc volt- 
age bias /i.2 are diagonal in energy with fi{£',£") = 
MS') 2n5{8' - f"), where = {e^£-^^)/T. + i)-i, 

Vi is the voltage applied, and is the electronic tem- 
perature. In the presence of time-dependent drive the 
occupation numbers fi{£',£") are not diagonal in en- 
ergy indices and do not commute. The voltage drive 
V{t) can be incorporated via the gauge transformation 
in time representation /i Ufv=oU^ with U{t',t") = 
e~*/o eV(t)dt^j-^/ _ ^//-j^ where we assume the convolution 
over internal time indices. 

Equation ^ can be simplified using polar decomposi- 
tion of the scattering matrijs^ 



S = 



U 
V 



S' 



U' 
v 



S' = 



R 
T 



(3) 

Here U, V, U', and V are unitary matrices in transport 
channels and T = diag(ri,T2, . . .) and R = 1 — T are 
diagonal matrices of transmission and reflection eigen- 
values. Equation ([T]) for the 2-terminal junction reduces 
to 



Six) = 2s Indet [l + f (S'A+S'A - l)] 

= 2 Indct + - ^) -/iA/TR(e'x _ r 

V/2VTR(e-*x _ 1) 1 + f^T{e'^ - 1) 



(4) 



Six) = Indet [1 + f (S+A^SA^ - l)] , 



(1) 



where 2^ takes into account the spin degree of freedom. 
Since the operators in the first (the second) row com- 
mute, the determinant can be calculated blockwise. Us- 
ing Eq. (j34p given in Appendix we obtain 
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Six) =2, ^ Tt£ ln[l + (1 - /i)/2Tp(e''^ _ i) + _ f^)Tp{e-'^ - 1)] 

=2, ^ Tr^ ln[l + h{\ ~ h)Tp{e'^ - 1) + (1 - f2)fiTp{e~'^ - 1)]. 
p 

I 



(5a) 
(5b) 



Here we used the matrix identity Indet(M) = 
Trln(M) and the trace Trg is taken in energy in- 
dices. The logarithm is taken assuming convolu- 
tion over internal energy indices, e.g., {fif2)s'S" = 

(27r)-i/rffi/i(f',fi)/2(fi,f")- 

In the following we show that the same result 
is obtained within the circuit theory of mesoscopic 
transport i^i^^i^^ In this case the system is represented 
by discrete circuit elements as depicted in Fig. [51 The 
cumulant generating function of the transferred charge is 
given byi^i^ 



J Tp / {Gi(x),G2(0)} _ J 
2 \ 2 



where 



G,(0) = 



1 2hi 
-1 



, (6) 



(7) 



are the quasiclassical Keldysh-Green's functions of the 
terminals with hi = 1 — 2/^. The counting field is incorpo- 
rated through the gauge transformation of the Keldysh- 
Green's function of the left terminal. 



Gi(x) =e-^>^^i/^ Gi(0) e'^^^i/^^ 



where 



1 

1 



(8) 



(9) 



in Keldysh space. The logarithm and the trace in Eq. ^ 
are taken assuming the convolution both in Keldysh and 
in energy indices. 

Equation ([6]) can be simplified by using matrix repre- 
sentation in Keldysh space. After rewriting Trln(- ■ ■ ) = 



C?i(0)\C?i(x) 

L / 



{Tp} 



62(0) 



FIG. 2; The schematic circuit-theory representation of a 2- 
terminal junction. The left and right terminals are character- 
ized by the quasiclassical Keldysh-Green's functions Gi,2(0). 
The counting field x is related to the charges which enter the 
left terminal through the cross-section indicated by dashed 
line. It can be incorporated into the corresponding Keldysh- 
Green's function via gauge transformation given by Eq. ((SJ. 



lndet(---) and taking the determinant by blocks using 
Eq. (l34l) . the result coincides with Eq. ([5]). This proves 
the equivalence of the circuit-theory expression for S{x) 
[Eq. ([5])] and the Levitov determinant formula [Eq. ([T|)] 
for a coherent 2-terminal scatterer. The circuit-theory 
approach becomes advantageous in the case of several 
scatterers in series^i^ or for a multiterminal mesoscopic 
conductor with large number of conduction channelsi^l 
It gives the charge transfer statistics in terms of the scat- 
tering properties of individual elements, thus circumvent- 
ing a non-trivial task of obtaining the scattering matrix 
of the composite system and averaging over the phase 
shifts. 

Equation ([5]) is valid for a dc bias applied and energy- 
dependent transmission probabilities. In this case the oc- 
cupation numbers /i_2 are diagonal in energy and com- 
mute with each other and with Tp[£). The trace over 
energy reduces simply to the integration and we obtain 

Six) = ^ 5] / rff In (1 + [1 - h[£)]h{£)Tp{e^'' - 1) 
^ p ■' 

+ h{£)[\- h{£)]Tp{e-^^-l)). (10) 

Here Iq is the measurement time which is the largest 
time scale in the problem, much larger than the charac- 
teristic time scale on which the current fluctuations are 
correlated.— The form of 5(x) reveals that the elemen- 
tary charge transfer processes are single-electron trans- 
fers to the left and to the right. The electron transfers at 
different energies and in different channels are indepen- 
dent. The term [I- fi{£)]f2{£)Tpe'^ in Eq. ^ describes 
the electron transfer from the right to the left lead at en- 
ergy £ in the channel p. The probability of this process 
is proportional to the probability /2(f) that the state in 
the right lead is occupied, the probability 1 — fi{£) that 
the state in the left lead is empty, and the probability 
Tp{£) of transfer across the scatterer. A similar analysis 
holds for the electron transfer from left to right. The left 
and right transfers are correlated because 5(x) is not a 
sum of the left- and right-transfer generating functions. 

Equation ([5]) is also valid for the time-dependent volt- 
age applied and energy-independent transmission proba- 
bilities. In this case fi^2{£' ,£") do not commute with 
each other. [However, the order of /i and /2 can be 
exchanged as shown by Eqs. ([5a|) and ((5b)) .] The loga- 
rithm has to be calculated with the matrix structure of 
/i^2 in energy indices taken into account because time- 
dependent drive mixes the electron states with different 
energies. 
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A few remarks on applicability of our approach are in 
place here. The circuit theory we have used applies to 
instantaneous scattering at the junction, with the fre- 
quency uj of the bias voltage much smaller than the in- 
verse dwell time r^"'^. It also applies for a dot sandwiched 
between 2 terminals and lo much smaller than the inverse 
RC time of the charge relaxation on the dot. In 
a typical experimental situation tuc T^r^^ with lo 
limited by the dwell time to <^ r^"^. However, as shown 
in Ref. [3^, the photon-assisted noise Sj{ijj) in the lead- 
ing order (eVo/w)^ <C 1 in driving amplitude Vq does 
not depend on t^: in a weakly -driven junction the fre- 
quency range is set by the RC time, oj <^ "^rc- 
dwell time becomes important for larger voltage ampli- 
tudes. For eVo/^ ^ 1 8' maximum in dSi/duj develops 
at the frequency ui ^ t^^ because the photon- assisted 
noise S'/(w) probes the electronic distribution function 
f{£;t) = Jdr f{t + T/2,t - t/2) e'^^ which relaxes on 
a timescale given by r^^^i^ At frequencies w ~ r^^^, 
the differential noise dSj/dVo as a function of Vq is in- 
creased with respect to the case lo <^ . The increase 
is of the order of 20% with the distances between maxima 
and minima only weakly affected. This is consistent with 
the experiments of Schoelkopf et al^ and Kozhevnikov 
et alS^ which can be described by assuming the instan- 
taneous scattering even though the frequency of the ac 
signal applied is comparable or larger than the inverse 
dwell time. 



III. ELEMENTARY CHARGE TRANSFER 
PROCESSES AT = 

A. Decomposition of the cumulant generating 
function into elementary processes 

In this section we study elementary charge transfer 
processes in a 2-terminal junction with constant trans- 
mission eigenvalues {Tp\ and time-dependent voltage 
V{t) applied. The cumulant generating function S{x) 
is given by Eq. ([5]) . To identify the elementary processes 
we diagonalize the operator under the logarithm in en- 
ergy indices. As a first step, we rewrite Eq. (|5bp in terms 
oi h= hi and h=h2'. 

5(x) = 2 ^ Trg ln[l -Tp{l-hh) sin2(x/2) 
p 

- iTpih - h) sin(x/2) cos(x/2)], (11) 

where h{£',£") = tanh(£:'/2Te)27r(5(£: ' - £") and h = 
UhU'^ ^ The S{x) further simplifies in the zero tempera- 
ture limit, in which the hermitian /i-operators are involu- 
tive h? = h? = I, and satisfy [/i, {h, h}] = [h, {h, h}] = 0. 
Thus the eigensubspaces of {h, h} are invariant with re- 
spect to h, h, and hh, and the diagonalization problem 
reduces to the subspaces of lower dimension. The typical 
eigensubspaces of {/i, K\ are 2-dimensional and spanned 



by the eigenvectors Vq and v_q, = hva of hh which cor- 
respond to the eigenvalues e^*" (a is real). The diago- 
nalization procedure in invariant subspaces is described 
in detail in Ref. '49'. 

For computational reasons it is convenient to impose 
periodic boundary conditions on the voltage drive V{t + 
t) = V{t) with the period t ~ 2tt/lo. In this case the 
operator h = UhU'^ couples only energies which differ by 
an integer multiple of lo. This allows to map the energy 
indices into the interval < £ < lo while retaining the 
discrete matrix structure in steps of lo. The operator hh 
in energy representation is given by 

{hh)nm{£) = {hh){£ + HLO, £ + muj) 

oo 

= sg\\{£ + nuj) ^ an+ka*ra+k sgn(f - kLO - eV), 

k— — oo 

(12) 

with 

a„ = - r dt e-'/o '^t'^AyctOg-^t. (13) 

T Jo 

Here V — {1/t) J V{t)dt is the dc voltage offset and 
^V{t) = V{t) — V is the ac voltage component. The 
coefficients a„ satisfy 

00 00 

an+ka*n+k = ^nm and ^ n|a„p = 0. (14) 

fc= — 00 n= — co 

The cumulant generating function S{x) at zero tem- 
perature is given by 5 = Sil + Sm + S2 with^ 

SilMx) = Ml,r 
X ^ ^ In [1 + TnRnPkLM {e'^ + e-'^ - 2)] , (15) 

n k 

and 

S2ix) = ^ E 1^ [1 + ^« (^"'''" - 1)] • (16) 

n 

Here Ml — toLOi/n, Mr — tQ{LO — uji) / it , and loi — eV — 
leV/uj\uj, where [xj is the largest integer less than or 
equal X. The coefficient k = ±1 in Eq. (fTB|) is related 
to direction of the charge transfer, with k = 1 (k = — 1) 
for eV > {eV < 0). The total measurement time is 
much larger than the period r and the characteristic time 
scale on which the current fluctuations are correlated. 

The parameters pkL,R in Eq. ([15]) depend on the de- 
tails of the time dependent voltage drive. They are given 
by PkL,R = sin^{akL,R/2), where e^'""^!") are the eigen- 
values of {hh)nm calculated for £ £ (0,a;i) [£ £ (0^1,0;)]. 
The eigenvalues are obtained using a finite-dimensional 
matrix {hh)nTm with the cutoff in indices n and m being 
much larger than the characteristic scale on which |a„| 
vanish. 
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V{t) = V + AV{t) 

FIG. 3: The elementary charge transfer processes in a voltage- 
driven junction. The unidirectional processes represent single 
electrons injected towards the contact due to excess dc off- 
set voltage V applied. The bidirectional processes represent 
electron-hole pairs created by time-dependent voltage drive 
and injected towards the contact. The probabilities pk of pair 
creations depend on the ac voltage component AV{t). 



Equations ([T5|) and give the charge transfer statis- 
tics in a 2-terminal junction driven by a time-dependent 
voltage. The result is valid in the zero-temperature 
limit, with no thermal excitations present. The prob- 
ability distribution of the number of charges N trans- 
ferred within measurement time is given by V{N) = 
(27r)~^/^^dx exp[5(x) — iiVx]- The cumulants of iV are 
given by C„ = [df^S{x)]x=o and are related to the higher- 
order current correlators at zero frequency ! "'^^'^^ For ex- 
ample, the average current, the current noise power, and 
the third cumulant of current fluctuations are given by 
/ = (e/io)a,x5|x=o, Si = (eVio)5fx'5|x=o, and C/ = 
(e3/to)9f-^5|x=o, respectively,^ 

The elementary charge transfer processes can be iden- 
tified from the form of S, similarly as in Refs.lsol andlsil. 
The result is depicted schematically in Fig.[3l The S^ix) 
describes unidirectional single-electron transfers due to 
the excess dc bias voltage V applied. The number of 
attempts for an electron to traverse the junction within 
measurement time to is given by to|el^|/7r (we take into 
account both spin orientations). The transfer events in 
different channels are independent. The term TnC"*"-^ in 
Eq. (|16p describes a single-electron transfer with proba- 
bility Tn in nth channel. The unidirectional charge trans- 
fer processes give contributions to the average current 
and higher-order cumulants. 

The Sil{x) and Sir{x) describe bidirectional charge 
transfer processes. Different bidirectional processes are 
labelled by k in Eq. These processes represent 

electron-hole pairs created in the source terminal by the 
time-dependent voltage drive and injected towards the 
scatterer. The probability of such an electron-hole pair 
creation is given by pk and depends on the details of 
the ac voltage component ISV{t). The charge transfer 
in nth channel occurs if one particle (e.g., electron) is 
transmitted and the hole is reflected, or vice versa. The 
probability for the whole process is given by T„RnPk- 



The bidirectional processes contribute to the noise and 
higher-order even cumulants. However, they give no con- 
tributions to the average current and higher-order odd 
cumulants because electrons and holes are transmitted 
with the same probability. 

The two types of bidirectional processes Sil and Sm 
differ in the number of attempts M^^n and have different 
probabilities PkL,R- The M^n depend on the number 
eV /uj of unidirectional attempts per period per spin. The 
simplest statistics is obtained for an integer value of eV /lo 
for which Sil vanishes, in agreement with Ref. 0. 

B. Comparison of different time-dependent 

voltages 

The elementary processes at zero temperature can be 
probed by noise measurements. In what follows we com- 
pare the elementary processes and the noise generated by 
different time-dependent voltages. We focus on standard 
periodic voltage signals such as cosine, square, triangle, 
and sawtooth. We also present results for Lorentzian 
voltage pulses which provide the simplest charge transfer 
statistics for certain amplitudes. The applied voltage is 
characterized by coefficients a„ given by Eq. (fT^ . These 
coefficients are calculated explicitly in Sec. B in Appendix 
for the driving voltages of interest. 

We first consider ac voltage drive with no dc bias ap- 
plied, V — Q. In this case only the bidirectional pro- 
cesses of R-type remain, S{x) = Sir{x)- The number 
of attempts during the measurement time is given by 
M = 2to/T which corresponds to a single attempt per 
voltage cycle per spin. The current noise power is given 
by 

with the probabilities pk of electron-hole pair creations 
obtained from the eigenvalues of {hh)nrm as discussed in 
SecHEl 

The probabilities pk for the harmonic voltage drive 
V{t) = Vocos{ujt) are shown in Fig. Ul^a). As the ampli- 
tude Vq increases the probabilities pk also increase and 
new pairs start to enter the transport. This results in the 
oscillatory change of the slope of Sj as a function of Vq . 
The decomposition of the differential noise dSj / dVo into 
elementary processes is shown in Fig. Iljb) . The differen- 
tial noise dSi/dVo for different time-dependent voltages 
is shown in Fig. [5] for comparison. 

The bidirectional processes with the unit probability 
Pk = 1 represent electron-hole pairs which are created and 
injected towards the scatterer in each voltage cycle. In 
this case the electron and hole transfers are statistically 
independent. This can be seen from the corresponding 
cumulant generating function which reduces to Snk = 
S^k + S^^k, where ^^.'^(x) = (2io/r) ln[l + T^{e^'>c _ i)] 
describe electron and hole transport. The electron-hole 
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FIG. 4: The probabilities pk of electron-hole pair creations 
for harmonic drive V{t) = Vo cos(ajt) as a function of the 
amplitude Vo [panel (a)]. With increasing Vo more and more 
pairs are created per voltage cycle. The decomposition of the 
differential noise dSi /d{eVo) (solid line) into contributions of 
elementary processes (dashed lines) is shown in panel (b). 



pairs which are created with the probability < < 1 
result in correlated electron and hole transfers. The 
corresponding cumulant generating function is given by 
Snkix) = (2io/T) ln[l+ T„i?„pfe (e*x + e-*x - 2)] and can- 
not be partitioned into independent electron and hole 
contributions. 

The interpretation of the shot noise in terms of 
electron-hole pair excitations has been studied previously 
in Ref. [H in the regime of low-amplitude harmonic driv- 
ing eVb/w <C 1. In this case only one electron-hole pair is 
excited per period with probability pi « (eVb/2w)^ ^ 1. 
Remarkably, Fig. H] shows that the single electron-hole 
pair is excited not only for small amplitudes but also for 
amplitudes comparable or even larger than the drive fre- 
quency. This extended range of validity can be covered 
by taking into account the higher-order terms in the ex- 
pression for probability pi « J2'^=i "-[^^nleVo/w)]^, where 
Jn{x) denote Bessel functions (cf. Sec. lIVp . The first 3 
terms approximate the exact pi shown in Fig. 3] to accu- 
racy better than 0.3% for eVo/w < 2. 

For large driving amplitudes the charge transfer statis- 
tics does not depend on the details of time-dependent 
drive^i^ and can be characterized by an effective volt- 
age Vcs — J_^dt\V{t)\ (here we assume V — 0). In 
this case there are N = |eKiff|/2a; ^ 1 processes with 
Pk — 1. The charge transfer statistics in the leading or- 
der in N consists of N^.h = 2N^/t uncorrelated elec- 
trons and holes injected towards the contact in each 
transport channel during measurement time: S{x) — 
iVe E„ ln[l + T^{e-'^ - 1)] + En + Tnie'"" - !)]■ 
This result can be interpreted by comparison with the 
generating function of a dc bias given by Eq. The 
electrons (holes) are injected during time intervals To 
(rh = r — Tc) per voltage cycle in which eV{t) > 
[eV{t) < 0]. For a large number of injected particles, the 
time-dependent drive in these intervals can be replaced 
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FIG. 5: The differential noise dSi /d{eVo) as a function of the 
amplitude Vo for different ac driving voltages {V = 0) is shown 
in top panel: square (dash-dotted), cosine (solid), triangle 
(dotted), sawtooth (dashed), and Lorentzian pulses (short- 
dotted line). The width of Lorentzian pulses is tl ~ O.lr. 
The oscillations are due to elementary processes which are 
created as the voltage amplitude increases (cf. Fig. |4ll . The 
decomposition of the differential noise into contributions of 
elementary processes is shown in lower panels. 



by effective dc voltages V^g and V^l^\ The number of 

attempts is given by Ni = {tQ/TT){n/T)\eV}^^\ {i = e, 
h), in agreement with Eq. The noise generated is 

the sum of independent electron and hole contributions 
Si = (le'^Vcff I/tt) ^„ T„i?„. This explains the asymp- 
totic behavior of differential noise at large driving ampli- 
tudes shown in Figs. [3] and O 

In the following we consider time-dependent voltage 
drive with a nonzero dc offset V ^ which creates 
both unidirectional and bidirectional elementary pro- 
cesses. The unidirectional processes generate dc noise 
which is given by Sdc = (e^/7r)(X]p Tp-Rp)et^. The bidi- 
rectional processes given by Eq. ([TS]) generate the excess 
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FIG. 6: The excess noise Si — Sdc is shown in panel (a) for dif- 
ferent time-dependent bias voltages; square (dashed), cosine 
(short-dotted), triangle (dash-dotted), sawtooth (dotted), 
and Lorentzian pulses (solid line). The width of Lorentzian 
pulses is TL ~ O.lr. The amplitude and the dc offset in- 
crease simultaneously, Vo = V. The probabilities pkL (dot- 
ted lines) and p^R (dashed lines) of bidirectional elementary 
processes are shown in panel (b) for the cosine voltage drive 
V{t) = V + V cos{LL}t). Only one elementary process of the L- 
type and one of the R-type can be created per period {k = 1). 
The decomposition of the excess noise Si — Sdc into contribu- 
tions of elementary processes is shown in panel (c). 



noise Sj - Sdc = (e'^ /ta)df Si\x=a: 



Si — Sd 



5Z "^p^p 




VPkL + (1 - V)P, 



'kR\ 



_(18) 

where v = eV /lo— [eV /lo\ is the fractional part of eV /u). 
For u 7^ 0, 1 there are 2 types of bidirectional processes 
(labelled by L and R) with different number of attempts 
and different probabilities of electron-hole pair creations. 

We study bidirectional processes first in the regime 
Vb = 1^ in which the ac amplitude and the dc offset 
increase simultaneously. The excess noise Si — Sdc for 
different time-dependent voltages is shown in Fig. [6l^a) 
for comparison. The probabihties PkL{R) ^ function of 
voltage are shown in Fig.[SI|b) for the cosine voltage drive. 
The decomposition of the excess noise into contributions 
of elementary processes is shown in Fig.lHlJc). For the co- 
sine voltage drive there are 2 bidirectional processes (one 
of L-type and another of R-type) which are excited per 



period. The L-type processes transform continuously into 
the R-type ones at integer values of eV/uj, while the R- 
type processes disappear. The step-like evolution of the 
probabilities of R-type processes as a function of voltage 
does not introduce discontinuities in the current noise 
because the corresponding number of attempts vanishes. 
Instead, the interplay between L- and R-processes results 
in kinks and the local minima at integer values of eV /lo, 
as shown in Fig. [6]^ a). The L-type (R-type) processes 
give the dominant contributions as eV /lo approaches the 
integer values from the left (right) because of the number 
of attempts which is proportional to v (1 — v). 



The Lorentzian voltage drive is special because it pro- 
vides the simplest one-particle charge transfer statistics. 
This is achieved for impulses carrying an integer number 
of charge quanta eVo/oo = N a.t offset voltages V >Vq. In 
this case the Lorentzian pulses do not create electron-hole 
pairs and only single-particle processes given by Eq. (fTB|) 
remain. The charge transfer statistics in each transport 
channel is exactly binomialM-'^ The noise is reduced to 
the minimal noise level 5'dc of the effective dc bias, as 
shown in Fig. E^a) . A formal reason for this is the van- 
ishing of the coefficients a„ = for n < —N [cf. Eqs. 
and (|42p ]. A many-body state which is created by opti- 
mal Lorentzian pulses has been obtained by Keeling et 
alM recently. 



Now we focus on regime in which the ac component 
AF(i) of the drive is fixed and the dc offset V changes. 
At low temperatures, the shot noise Si is a piecewise 
linear function of the dc voltage offset V with kinks cor- 
responding to integer multiples of the driving frequency 
eV/Lo = N and slopes which depend on the shape and 
the amplitude of the ac voltage component i^i^i The dif- 
ferential noise dSi/d{eV) consists of a series of steps, 
as shown in Fig. [7]^ a). The piecewise linear dependence 
of noise can be understood in terms of probabilities of 
elementary processes.^ For the ac component fixed, the 
probabilities pkR and pkL of electron-hole pair creations 
are piecewise constant as a function of V , and can be re- 
labelled by pkB.{y) = Pfc^^ and PkhiV) = Pfe^'''^' where 
N = \_eV /uj\ . This results in the piecewise linear depen- 
dence of the excess noise given by Eq. (fT5|) as a function 
of V. For the dc offset in the interval m — 1 < eV /lo < m 
with m integer, the excess noise is a linear combination of 
and p^,™'' (fc = 1, 2, . . .) with the number 
V and V, respectively. For 



processes pf, 
of attempts proportional to 1 



and 



TO < eV /lo < to 4- 1 a different set of processes p^" 

pI™^^'' contributes and the excess noise changes slope. 
The excess noise for harmonic drive is shown in Fig. [7jb) 
(solid line) with decomposition into contributions of el- 
ementary processes (dotted lines). The probabilities of 
elementary processes are shown in the inset. 
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FIG. 7: The differential noise dSi / d{eV) as a function of dc 
offset V is shown in panel (a) for different time-dependent bias 
voltages: cosine (solid) , square (dashed) , and sawtooth (dash- 
dotted line). The ac amplitude of the drive is kept constant, 
eVo /io = ^. The excess noise Si — Sdc as a function of V (solid 
line) and the decomposition into contributions of elementary 
processes (dotted lines) are shown for the cosine drive in panel 
(b) [cf. Eq. (|18|) ]. The probabilities of elementary processes 
are shown in the inset. 



C. Beam splitter geometry 

Here we study a multiterminal beam splitter geometry 
depicted in Fig. [S] The source terminal is biased with 
a time-dependent periodic voltage V{t) and through a 
mesoscopic conductor attached to several outgoing ter- 
minals. The conductor is characterized by a set of trans- 
mission eigenvalues {Tp} and the outgoing leads by con- 
ductances gi. We are interested in the limit in which 
the outgoing leads play a role of a detector and only 
weakly perturb the charge transfer across the conduc- 
tor. This is achieved for the conductance g^: = J2i 9i 
to the outgoing leads much larger than the conductance 
g = (e^/7r)^pTp of a conductor. In this case the par- 
ticles which traverse the conductor enter the outgoing 
terminals with negligible backreflection into the source 
terminal. A similar setup with spin-selective outgoing 
contacts has been used in Ref. [6l| to reveal singlet elec- 
tron states. 

The cumulant generating function is calculated using 
the circuit theory similarly as in Sec.|TT]for the 2-terminal 
case. In contrast to Sec. [ITl here we assign the counting 




GiiXi). 



FIG. 8: Schematic representation of a multiterminal beam 
splitter. The source terminal is biased with a time-dependent 
voltage V{t) with respect to the outgoing terminals. The total 
conductance gj: — gi of the outgoing leads is assumed to 
be much larger than the conductance g = (e'^/yr) Tp of the 
source contact. 



fields Xi to the outgoing terminals. The Green's function 
G{Q) of the source terminal and the Green's functions 
GiiXi) of the outgoing ones are given by 



G(0) = 



1 2h 

-1 



and 



2h 

-1 



(19a) 



(19b) 



Here h and h are the matrices in energy indices defined in 
Sec. nil Al The Green's function Gc of the internal node is 
given by matrix current conservation and normalization 
condition = 1. In the limit ^ g, the central node 
is strongly coupled to the outgoing terminals and the Gc 
can be obtained in the lowest order with the terminal G 
unattached. For simplicity we assume tunnel couplings 
to the outgoing terminals. In this case matrix current 
conservation reduces to—i^ 



^giGi{Xi),Gc 



= 0. 



(20) 



In the following we work in the low-temperature limit 
in which = ~ 1. We seek for the solution in the 
form Gc = Pc T where f = (fi,f2,'f3) is the vector 
of Pauli matrices in Keldysh space. From Eqs. (fTO)) 
and (|20p we obtain Pc = {h,ihc — s,ihs + c) where 
c = J2t9iCOs{xz), s = I]i5»sin(xi), and = gi/g^. The 
cumulant generating function of the charge transferred is 
given by 



Tr In 



Tp / {G(0),Gc({x})} 
2 V 2 



(21) 

where summation over the internal Keldysh and energy 
indices is assumed. The diagonalization of {G, Gc}^,^,, 
is performed along the lines of Ref. H^. 
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For the cumulant generating function we obtain 
iS({x}) — SiL + Sin + '^2 where Si^^n are the contribu- 
tions of bidirectional processes and 1S2 is the contribution 
of the unidirectional ones. The unidirectional processes 
are described by 

S2{{x}) = ^E^'^ (^1 + T„E.9Ke''^' - 1)) - (22) 

where to is the measurement time and V is the dc volt- 
age offset (here we assume eV > 0). The bidirectional 
processes are described by 



xEEin 



+Pk 




PkaTnRn 



2) 



5*5j(e 



2) 



(23) 



Here k labels the bidirectional processes, n labels trans- 
port channels, i and j label the outgoing terminals, and 
a = L, R. The number of attempts and the proba- 
bilities pka are the same as in Sec. IIII Al 

Equations and give the statistics of the charge 
transfer in a multiterminal beam splitter at low temper- 
atures in the presence of a periodic time-dependent drive 
at the source terminal. The Si{{x}) and S2{{x}) have 
a direct physical interpretation. The unidirectional pro- 
cesses, which are described by 52, are the single-electron 
transfers across the structure due to the dc offset voltage 
V applied to the source terminal. The term r„gje'-^' in 
Eq. (j22p represents the process in which an electron in 
the nth transport channel traverses the conductor with 
probability T„ and enters the outgoing terminal i with 
probability gt = gi/gi:- 

On the other hand, the bidirectional processes repre- 
sent the electron-hole pairs created in the source terminal 
and injected towards the conductor. The probabilities of 
such excitations are given by pka ■ The interpretation of 
the cumulant generating function given by Eq. (j23p can 
be obtained from a simple counting argument. The term 
PkTnRngiS^-^' in Eq. (|23p represents the process in which 
electron-hole excitation is created, hole is reflected, and 
electron is transmitted into the outgoing terminal i. Sim- 
ilarly, the term p/cT„i?„gie~*-^' represents the process in 
which the electron is reflected and the hole is transmitted. 
Finally, the term PkT^gigje^^' e~'^^^ {i 7^ j) represents 
the process in which both particles are transmitted with 
electron entering terminal i and hole entering terminal j . 

Charge transfer statistics and current correlations can 
be obtained using the cumulant generating function 
S{{x})- For example, the current cross correlation be- 
tween different terminals i and j is given by Sij = 
(eVto)C,.x.'5|x=o: 



Si' 



E^n 9^~9JJ2^VpkL + {l-v)pkR]■ (24) 



We find that only bidirectional processes in which both 
particles are transferred (one into the terminal i and an- 
other into the terminal j) give the contributions to the 
cross correlation 5*^. This has been obtained previously 
by Rychkov et alfi- in the special case of harmonic drive 
in the limit of small driving amplitudes (V" = and 
eVo/uj < 1). 

The cross correlation Sij depends on bidirectional pro- 
cesses and is proportional to the excess noise Sj — Sdc of 
a 2-terminal junction [Eq. (jl8[) ]. In a 2-termianal setup, 
the excess noise is just a small correction to Sdc gener- 
ated by unidirectional processes for a bias voltage with 
y ~ Vb (Fig. [5]). The contribution of Sdc component can 
be reduced by measuring current cross correlations be- 
tween different outgoing terminals in the beam splitter 
geometry with negligible backscattering. 



IV. CUMULANTS AT Te 7^ 

The full counting statistics and the corresponding el- 
ementary transport processes obtained in Sec. IIIII are 
valid description in the low temperature limit only, in 
which electron-hole pairs are created by the applied volt- 
age and no thermally excited pairs exist. Formally, the 
diagonalization of the operator in Eq. ([TT|) in energy in- 
dices, which is needed to deduce the elementary pro- 
cesses, is based on the involution property of /i-operators: 
h? = h? = 1. This property no longer holds at finite tem- 
peratures which are comparable to the applied voltage. 
Nevertheless, the method we use enables the efficient and 
systematic analytic calculation of the higher-order cumu- 
lants at finite temperatures. 

The cumulants can be obtained directly from Eq. (fTT|) 
by expansion in the counting field to the certain order 
before taking the trace. The trace of a finite number of 
terms can be taken in the original basis in which h and h 
are defined. In the following we illustrate the approach 
by calculation of the average current, the current noise 
power, and the third cumulant at finite temperatures. 
From Eq. pT|) we obtain 

9,x5|x=o-E^p'^^'£(^"'^)' (25a) 
p 

dj^S\^=o ^Y.{Tj,Trs{l - hh) - (TpV2) TTe[{h - hf]}, 



(25b) 



9f^5|;,=o=E{(^p/2)Tr£(/i=^-/»') 



+ (3/2)rp2(i _ Tp) Tre[hh{h - h)] 

+ rp[l- (3V2)]Tr£(/i-/i)}. (25c) 

In energy representation, the operators h and h are given 
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by 

h{£',£") ^h{£') 2TrS{£'-£") 

= akal^^h{£') 2nS{£" - £' - mu)), (26a 



h{£', £") = ^kal+,M£' -kto- eV) 

k.7n 

X 2ttS{£" -£' ~muj). 



(26b) 



Here h{£) — tanh(f/2re) and we used the properties 
of {a„} given by Eq. (HH). After integration over en- 
ergy in Eqs. we obtain the average current / — 
(e^/7r)( Tp)V. The current noise power and the third 
cumulant are given by 



Si 



p p 

eV + nuj 



\anf {eV + nuj) coth 



n— — oo 



2n 



and 



^ p 

E 



CXD 

X > |a,; 

n— — OO 



2rp coth 



+ nuj 



2T, 



'{eV + nuj) coth 



2 / eV^ + nuj 



2T, 



(28) 



respectively. For a dc voltage bias a„ = (5„.o and the 
noise and the third cumulant reduce to 5/ — eFI and 
Ci = e^F^I at zero temperature. Here F — ~ 
Tp)]/iEpTp) is the Fano factor and F3 = [EpTp{l - 

Tp){l-2Tp)]/{Y.pTp). 

The average current is linear in dc voltage offset, 
which is consistent with the initial assumption of energy- 
independent transmission eigenvalues and instant scat- 
tering at the contact. The result for the current noise 
power, Eq. (|27p . describes the photon-assisted noise for 
arbitrary periodic voltage drive. The coefficients a„ for 
harmonic drive V{t) — V + Vo co&{ujt) are given by the 
Bessel functions, a„ — JnieV^lLS), and Eq. (P7)) reduces 
to the previous results obtained by Lesovik and LevitOA*^ 
and Pedersen and Biittiker— (see also Ref. [ll). The 
accurate noise measurements at finite temperatures in 
the presence of the harmonic driving are performed in 
Ref. [S^. The results are in agreement with Eq. ((27|) . 

In the following we discuss the low- and high- 
temperature limits of Si and C/. At high temperatures 
Te ^ |ey I, rtow, with n^uj being the characteristic energy 
scale on which |a„g| vanish, the current noise power re- 
duces to the thermal equilibrium value Si = 2TeG, which 
is just a manifestation of the fluctuation-dissipation the- 
orem. The third cumulant is in this regime proportional 



to the average current, C/ = e^FI. At high tempera- 
tures Si and Ci carry no information on the details of 
the time-dependent voltage drive. 

At low temperatures ^ uquj, the current noise 
power reduces to 



Si = -iY,TpRp) J2 \eV + nij\\a,,\\ (29) 



The differential noise dSi/dV is a piecewise constant 
function of V with steps given by^ 



A 



dSi 
W 



eV/ui=N 




N-0 



(30) 



All steps add up to {2e^/Tr)^ TpRp because of 



(27) EJ 



= 1, see Fig. ma) ,52. The steps in dSi/dV 



have been measured for harmonic drive in normal^^ and 
normal-superconductor^ junctions. In the supercon- 
ducting state they appear at integer values of 2eV/ui, 
which can be interpreted as a signature of the elemen- 
tary charge transport processes in units of e* = 2e. 
The effective charge is doubled in the superconducting 
state due to the Andreev processes. We point out that 
for a general voltage drive, certain steps at integer val- 
ues of e*V/uj = N may vanish if the corresponding co- 
efficient a_Ar = 0. For example, for a square-shaped 
drive with integer amplitude e*Vo/uj = m, the steps 
a.t e*V/uj ^ m + 2k (k ^ 0) vanish [cf. Eq. §7} and 
Fig.Efa)]. 

At low temperatures, the third cumulant reduces to 
Ci = e^F^I. Unlike the current noise power, the third 
cumulant at low temperatures does not depend on the 
ac component of the voltage drive. This is because the 
bidirectional processes, created by the ac voltage compo- 
nent, do not contribute to odd-order cumulants at low 
temperatures [recall Eq. (fTSjl]. 

We conclude this section by comparison of two formu- 
las for the current noise power at zero temperature. For 
simplicity we take V = 0. Equation ([M)) for the current 
noise power reduces to 



Si 



2e^uj 



(31) 



Here we used Eq. to restrict the summation to the 
positive n only. On the other hand, the current noise 
power is also given by Eq. pT|) . Regardless the simi- 
lar form, the physical content of these two equations is 
very different. Both equations give the same result for 
Si as a consequence of the invariance of trace. However, 
Eq. (fT7|) has been obtained by taking the trace of the 
operator in Eq. (|lip in the basis in which it is diagonal. 
The cumulant generating function given by Eq. (jlSp is 
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decomposed into contributions of elementary and statis- 
tically independent processes. The terms proportional to 
Pk which appear in Eq. (fTT]) are the contributions of these 
processes to the noise. Equation (|3ip has been obtained 
by taking the trace of the operator in Eq. (fTTj) in a ba- 
sis in which it is not diagonal. Although the end result 
for Si is the same, the individual terms proportional to 
n|a„p which appear in Eq. (|3ip have no direct physical 
interpretation. 

In the limit of small-amplitude voltage drive, eVo/w 
1, only one electron-hole pair is excited per period with 
probability pi. Comparing Eqs. pT|) and (Plj) we obtain 
that in this case pi = X^^i'^l'^nP- ia-ct, as shown 
in Fig. |4] for harmonic drive, the assumption of small 
amplitudes can be relaxed to the amplitudes comparable 
or even larger than the drive frequency. The accuracy of 
this approximation depends on the actual voltage drive 
(cf. Figs. H and E]). 



V. CHARGE TRANSFER PROCESSES AT / 

Here we study the effect of finite temperature on ele- 
mentary charge transfer processes. For a periodic voltage 
drive, the cumulant generating function S{x) in Eq. (jlip 
can be recast in the following form: 



0.25 



0.25 



to 



E 



dE In 1 



-1) 



(32) 

Here qnG^^^ describes the process in which n charges are 
transmitted through the scatterer with probability g„. 
The charge transfer can occur in either direction depend- 
ing on the sign of n. To obtain the probabilities (/„ we 
diagonalize the operator under the logarithm in Eq. (Ilip 
numerically, multiply all the eigenvalues, and take the in- 
verse Fourier transformation with respect to x- For sim- 
plicity we focus on ac voltage drive with no dc offset. For 
a given ac drive, the probabilities (/„ depend on energy, 
temperature, and transmission, qn = qn{£ /uj,Te/uj;T). 
Transport properties of the junction enter through sum- 
mation over the transmission eigenchannels {Tp}. 

Equation ((32|) allows us to study the crossover from ac- 
driven fluctuations at zero temperature to thermal fluc- 
tuations at temperatures much larger than the voltage 
drive. The important difference between the two regimes 
is that ac drive mixes electron states of different ener- 
gies in contrast to thermal fluctuations which are diag- 
onal in energy. As we have seen, this is also reflected 
in the charge transfer statistics. The statistics for an ac 
drive (at zero temperature) can be interpreted in terms 
of electron-hole pairs after mapping the problem into en- 
ergy interval set by the driving frequency lu. The mixing 
of energy states is taken into account by diagonalization 
of the remaining matrix structure in energy. On the other 
hand, the statistics in the thermal limit is composed of 
electron transfers (in either direction) which are indepen- 
dent at different energies. In the crossover regime there is 
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FIG. 9: The probabilities q±n of multiple charge transfers 
for harmonic voltage drive V{t) = Vb cos(tjt) and transmis- 
sion Tp = 0.5. The temperature dependence of q„ is shown 
in panel (a) for different driving amplitudes. The energy de- 
pendence of q„ at low temperatures is indicated by plots for 
£/uj = (top line) and £/uj = 0.5 (bottom line). The energy 
dependence of q„ is shown in the thermal limit eVo/co = 
(b) and for eVo/co = 1 (c). At zero temperature q±i do not 
depend on energy and are given by q±i = T{1 — T)pi{eVo/Lu) 
(cf. Fig. As the temperature increases, multiple charge 
transfers start to enter the transport. The limit of thermal 
fluctuations is reached at temperatures larger than the driving 
amplitude. It can be recast again into single-charge transfers 
independent in the full energy range. 



still some mixing of different energy states present and we 
use the probabilities qn in Eq. ([321) with energy mapped 
in the w-interval. These different physical situations are 
depicted schematically in Fig. [TJ 

The probabilities (?„ are shown in Fig. [9] for cosine 
voltage drive V{t) — Vq cos{u}t) and a transport chan- 
nel of transmission Tp — 0.5. We first discuss thermal 
limit of probabilities <?„ shown by solid lines in Figs. [5] 
(a), (b). For w 3> Tg, the mapping into interval lo is 
irrelevant since it is larger than the energy scale of ther- 
mal fluctuations. In this case the probabilities reduce to 
q±i — /(I — f)Tp and g„ = otherwise, in accordance 
withEq. ifTO)) [cf. Fig.[9](b)]. As the mapping interval be- 
comes comparable or smaller than the temperature, the 
higher-order "bands" qn {n = ±1,±2,...) start to ap- 
pear. For uj < 2Te the probabilities qn are independent 
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of energy [Fig. [9] (a)] and are given by 

g„ ^ r ^ e^^^)-^-^ (33) 

with S{x) = -4(Te/w) arcsin2[y7^sin(x/2)]. Here S{x) 
is related to the cumulant generating function S{x) in 
thermal equilibrium by S{x) = (S^o/t) X)p'^(x)- The 
latter is obtained from Eq. (fTO|) after energy integration;^ 

We emphasize again that the picture of multiple charge 
transfers in thermal equilibrium limit is solely due to 
mapping into the energy interval u) in Eq. ([5^ . Although 
such mapping is not needed because thermal fluctuations 
are diagonal in energy, we nevertheless perform it here 
to provide a limit for ac-driven processes as the temper- 
ature increases. At temperatures much larger than ac 
drive there is no mixing of different energy states and 
the statistics in the full energy range reduces to single 
electron transfers which are independent at different en- 
ergies [Eq. pO))]. 

Now we focus on the ac-driven processes and the 
crossover region between ac and thermal fluctuations as 
the temperature increases. For simplicity, we consider 
driving amplitudes eVo/w = 1 and eVQ/uj = 2. In this 
case q±i{Tp) — Tp{l ~ Tp)pi(eVQ / uj) at zero temperature 
and do not depend on energy. These single-charge trans- 
fers originate from an electron-hole pair which is created 
with probability pi per voltage cycle (Fig. 2]). As the 
temperature increases, the interplay of thermal and ac 
excitations introduces a nontrivial energy dependence of 
probabilities q±i shown in Fig. [5](c). It also modifies the 
dependence on transmission eigenvalues with q±i{Tp) no 
longer being proportional to Tp(l — Tp) [cf. Eq. l|33|) ]. 
As the temperature increases further, the multiple-charge 
transfers come into play as shown in Fig. [5] (a) . At tem- 
peratures larger than the driving amplitude the limit of 
thermal transport is reached with statistics which can be 
recast again in the form of single-charge transfers, inde- 
pendent at different energies. 



correlators. The bidirectional processes contribute to the 
noise and higher-order even cumulants of the transferred 
charge at low temperature and give no contributions to 
the average current and higher-order odd cumulants. For 
an ac voltage with no dc offset, the noise is entirely due 
to bidirectional processes. The individual processes can 
be identified from the oscillations of the differential noise 
dSi/dVo as the amplitude Vq of the drive increases. 

A time-dependent voltage drive with a nonzero dc off- 
set generates both unidirectional and bidirectional pro- 
cesses. The bidirectional processes give rise to the excess 
noise with respect to dc noise level which is set by the 
unidirectional processes. The excess noise can be probed 
by measuring current cross correlations between different 
outgoing terminals in the beam splitter geometry. The 
cross correlations are only due to processes in which the 
incoming electron-hole pair is split and the particles enter 
different outgoing terminals. 

The method we use enables the systematic calculation 
of the higher-order current correlators at finite tempera- 
tures by expansion of the cumulant generating function 
in the counting field. We have obtained the current noise 
power and the third cumulant for arbitrary periodic volt- 
age applied. 

We have also studied the effect of finite temperature 
on elementary charge transfer processes. In the limits of 
low (high) temperatures the single-charge transfers occur 
in either direction due to ac excited (thermally excited) 
electron-hole pairs. However, the nature of the two is 
very different: the ac drive mixes the electron states of 
different energies while thermal fluctuations are diagonal 
in energy. In the crossover region there is still some mix- 
ing of different energy states present. Such an interplay 
of thermal and ac excitations can be interpreted in terms 
of multiple-charge transfers with probabilities dependent 
on energy and temperature. 
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APPENDIX 

A. Determinants of block matrices 

Let A, B, C, and D be the quadratic matrices of the 
same size. Then the following equalities hold: 



det 



A B 
C D 



det (AD - 


BC), 


[C,D] 


= 0, 


det (DA - 


BC), 


[B,D] 


= 0, 


det (DA - 


CB), 


[A,B] 


= 0, 


det (AD - 


CB), 


[A,C] 


= 0. 



(34) 
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In the case in which more than two blocks commute with 
each other, the corresponding determinants on the right 
hand side of Eq. (p4| coincide. 



The triangle voltage drive is characterized by ^V{t) = 
Wot/r - t/o for < i < t/2 and AV{t) = -Wot/r + 3Vb 
for t/2 < < < r. In this case 



B. Coefficients a„ 

Here we calculate coefficients {a„} given by Eq. (fT3|) for 
different time-dependent voltages. For the cosine voltage 
drive V{t) = V + Vocos{ujt), the coefficients a„ can be 
calculated using Jacobi- Anger expansion>^ 



^iz sin(^) 



(35) 



n— — oo 



where J„ are the Bessel functions of the first kind. From 
Eq. p3|) we obtain 



a-n = JnieVo/uj). 



(36) 



The square voltage drive is given by AV{t) — Vq for 
< t < t/2 and AV{t) = -Vo for t/2 < t < t. For 
noninteger values of eVo/oj, the coefficients a„ are given 
by 

^ 2 eFo sin[(n- eVb/^)7r/2] ^„_eVo/u)n/2 (nr\ 
TT w rt^ — (eVo/w) 

For integer values of eVb/i^, the coefficients a„ are ob- 
tained by taking the limit in the previous formula. 

The sawtooth voltage drive is given by AV{t) — 
2Vot/T - Vo for < t < r. In this case 



1 / TT (eVo/uJ + ~ eVo/2uj 

— , exp I ; 

2^2eVo/u} \ 2 eVo/w 



erf 



V^(eVb/a;- n) ,^/4 



erf 



V2eV^ / 



(38) 



where erf(z) = (2/^) dt exp(— i^) is the error func- 
tion. 



1 ^ r f.7T{eVo/uj + n)^-eVo/u; 

— , Ke < exp I — --— 

2yiTV^ I V 4 eVo/uj 



erf 



^/^{eVQ/uj~n) ^,^/A 
2y/eW^ J 

erf ( )f^^^^hi±A,^-!^ 



, . , , - (39) 

2^fW^ ' ' 

The voltage drive which consists of Lorentzian voltage 
pulses of width is given by 



AV{t) 



V- TTL 

TT ^ (t~ kTf + t] 



Vo sinh(27rTL/r) 



cosh(27rrL/r) — cos(27ri/r) 



(40) 



The total voltage V{t) — V + AV{t) varies between 
^min = V + Vb[tanh(7rTL/r) - 1] and Kiax = V + 
Vb [cotanh(7rri /r) — 1] . The coefficients a„ are given by 



, {sin[7r(a; + ^y)]}^^«/'' 
_i/2 {sin[7r(a; - iy)]Y^o/<^ 

y i2-n(eVo/u}+n)x 



, (41) 



where y — tj^/t. For integer values eVo/w = -/V > we 
obtain the simplified expressions: 



(-l)^27ri Res 



sin \ii{x + iy)\ i27r(n+w)a 



^=iv Vsin [7r(x — iy)] 

(42) 

for n > an ~ for n < —N, and a„=_jv = 

(—1)^6^^^^^^. In particular, the coefficients a„ for 
eVo/uj = 1 are given by a„ = g-^irny _ g-27i-(n+2)i/ 

— 1, a_i = — e^^'^^, and a„ = otherwise. For eVo/ui = 2 
we obtain a„ = e-2™a(i _ e-4'^2')[n +l-{n + 3)e~'^''y] 
for n > —2, a_2 = e^"*^^, and a„ = otherwise. 
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